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If the positron fraction and combined electron-positron flux excesses recently observed by 
PAMELA, Fermi and HESS have a dark matter origin, final state radiation (FSR) photons from 
dark matter annihilation into lepton-rich final states may be detected with observations of satellite 
dwarf galaxies of the Milky Way by ground-based atmospheric Cherenkov telescopes (ACTs). 
We find that current and near-future ACTs have excellent potential for such detection, although a 
discovery cannot be guaranteed due to large uncertainties in the distribution of dark matter within 
the dwarfs. We find that models predicting dark matter annihilation into two-Iepton final states 
and those favoring four-Iepton final states (as in, for example, "axion portal" models) can be re- 
liably distinguished using the FSR photon spectrum once measured, and the dark matter particle 
mass can also be accurately determined. 
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1. Introduction : Final State Radiation, Dwarf Galaxies, ACTs 



Recent measurements of the positron fraction in cosmic rays in the 10-80 GeV range by 
PAMELA ^ and the combined electron-positron flux up to the TeV scale by Fermi [Q] and HESS 
[^] indicate excesses inconsistent with conventional astrophysical background. These excesses can 
be explained, among other possibilities, by dark matter with an annihilation cross section of (av) ~ 
3 X 10~^^ cm^s~^ in the Milky Way and a mass of 1-3 TeV, provided annihilation is predominantly 
into electrons or muons |Q ^. 

If these excesses indeed have a dark matter origin, accompanying signals are expected in the 
form of energetic gamma rays. For leptophilic dark matter, the high energy end of the gamma 
ray signal is dominated by final state radiation (FSR). Dwarf galaxies — made up almost entirely 
of dark matter, with no detected neutral or ionized gas, minimal dust, no magnetic fields, and 
little or no recent star formation activity — are favorable targets for searches for such gamma 
rays. And since dark matter gamma ray signals from dwarf galaxies are mainly constrained by 
low statistics, atmospheric Cherenkov telescopes (ACTs), with typical effective areas ~ 10^ times 
that of Fermi, offer a distinct advantage. In this work, we focus on the prospects of detecting FSR 
from dark matter annihilation with ACT observations of dwarf galaxies. We also investigate the 
important question of whether the ACTs can measure the FSR photon spectrum precisely enough 
to distinguish between different leptophilic dark matter models. For details, please refer to [^. 

Final state radiation (FSR) — FSR is present whenever dark matter annihilates to charged 
particles, as is the case for the following three leptophilic models motivated by fits to PAMELA, 
Fermi, and HESS data ^ (with dark matter denoted by ;|^): 
(i) Model A: — )- (ii) Model B: 4e. (iii) Model C: 

^ denotes an intermediate "portal" particle, with mass taken to be of order 1 GeV. Factorization the- 
orems ensure that the energy spectrum of the FSR photons is, to leading order in m/ /m^^, indepen- 
dent of the details of the annihilation process, allowing for quasi-model-independent predictions. 
For annihilation into a lepton-antilepton pair II as in model A, the FSR flux within the leading log 
approximation is [0] 
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(1.1) 

(1.2) 



Here x = lEyj ^/s = Ey/ni),, <I>o = 1.4 x lO^''^ cm^^s^'GeV and 7 is a dimensionless factor 
that carries all the astrophysics information;^ astrophysical uncertainties, therefore, do not affect 
the spectrum of final state radiation. The spectrum in Eq. dEl] ) features a characteristic "edge" at 
the dark matter mass |^. 



' Eq. (1.1) and similar formulas below are equally applicable for decaying dark matter, with an appropriate redef- 
inition of the J factor. However, we do not consider decaying dark matter in this paper since the resulting FSR signals 
from dwarf galaxies are generally too weak to be detected. 
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For 4-body annihilation as in models B and C, the FSR spectrum is 

^ (i^) (100GeVy^2-x + 2xlogx-x^ / ^\ ^ ^^ 3^ 

dx \\pb J \ J X ) 

For models A and C, annihilation is to muons, and final state radiation form the subsequent decay 
of the muon should also be taken into account. Relevant formulas for this contribution are as 
given in [||, Typically, FSR off muons from the annihilation process remains dominant unless 
~ m^[|8|]. In this work we fix our parameters to (av) = 3 x 10^^^ cm-'s^^ m^^ = 3 TeV, and 
= \ GeV, as favored by fits to PAMELA, Fermi, and HESS data. It should be kept in mind 
that lower velocities in dwarf galaxies can lead to larger cross sections via greater Sommerfeld 
enhancement. 

Dwarf galaxies — We use the following dwarf galaxies, which have been known to be promis- 
ing candidates for dark matter searches, in our analysis: 

Draco (18.63 ± 0.60) Ursa Minor (18.79 ± 1.76) 

Willman 1 (19.55 ± 0.98) Segue 1 (19.0 ± 0.6) 
The number in parenthesis is logjo(Lx GeV^^cm^), where L is the astrophysical factor as defined 
in Eq. (L2) and calculated in [|8[ ^ ^; it should be noted that the uncertainties on these astrophys- 
ical factors are extremely large at present. The Sloan Digital Sky Survey has recently discovered 
many new dwarf galaxies, and since only a small region of the galactic neighborhood has been 
completely surveyed, several hundred more low-luminosity, dark matter dominated dwarf galaxies 
might still be discovered in the future. It would be straightforward to apply the analysis of this 
paper to any promising new dwarf that may be discovered, once the distribution of dark matter is 
mapped out to allow for at least an approximate determination of its L factor. 

Atmospheric Cherenkov telescopes (ACTs) — The key parameters governing the sensitivity 
of an ACT in observations of dwarf galaxies are its effective area A^ff, energy resolution e, and 
energy threshold. There are several currently operational (eg. MAGIC, VERITAS) and near future 
(eg. CTA) ACTs relevant for indirect dark matter detection; we refer to ||^ for more details, 
references, and individual key parameter values. Typically, current ACTs have Agff ~ 10^ cm^ 
and £ ~ 0.15, while future instruments are expected to reach A^ff ~ lO'^cm^ and e ~ 0.10. The 
typical instrumental energy threshold for ACT's is about 200 GeV. 

Previous observations and upper bounds — The dwarf galaxies mentioned above have been 
observed by ACTs without any positive detection, resulting in upper bounds on high energy gamma 
ray flux from dark matter annihilation or decay in these galaxies. Figure H, which compares these 
experimental bounds with theoretical predictions from the three leptophilic models of interest in 
this paper, shows that the predictions are consistent with the observed null results within the uncer- 
tainties in the astrophysical factors; for more details on bounds from individual observations, the 
reader is referred to ^ and references therein. 

Backgrounds — Since the dwarf galaxies themselves are not expected to contain significant 
sources of hard gamma rays of astrophysical origin, the background can be effectively measured 
by looking at a region of the sky close to the dwarf (called the OFF region); subtracting the OFF 



^The astrophysical factor for Segue 1 listed here represents an updated value |g] that was not available at the time 
of writing of [^]; all plots and discussions in the following sections use this updated value. 
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Figure 1: Comparison of experimental bounds with predictions from theory. The horizontal lines repre- 
sent the bounds from experimental observations ([^ and references therein). The three vertical bars for 
each search are the corresponding predictions of models A (left bar), B (center), and C (right), using the 
astrophysical factors as listed in the text, with circles denoting central values. 

region flux from the flux in the ON region (which contains the dwarf) eliminates the background up 
to statistical fluctuations. Astrophysical backgrounds can be estimated with standard extrapolations 
of charged lepton, hadron, and gamma ray spectra [ jl0| , 0] ; the hadronic background depends on 
the hadron rejection capabilities of the instrument ( see ^ for a more detailed treatment). Gamma 
rays from dark matter annihilation within the Milky Way also contribute to this background. In 
addition to FSR, these come from inverse Compton scattering (ICS) of starlight and CMB photons 
off energetic leptons from dark matter annihilation; we estimate the ICS contribution using the 



semi-analytic formalism in [|11[]. 

2. Detection Prospects 

Standard requirements for detection are (i) the excess in the ON region relative to the OFF 
region exceeds 3a (or 5a), and (ii) more than 25 (signal) events (after background subtraction): 

Significance = ^^'^''^^ _ . > 3 (or 5), (2.1) 

Number of signal events = OyAgfff > 25, (2.2) 

Here t is the observation time, denotes gamma ray flux, and Q. denotes solid angle. For our esti- 
mates we ignore systematic errors and assume background subtraction with ON and OFF regions 
of the same size. To improve sensitivity, it is useful to choose an energy threshold that maximizes 
the ratio <I>signai/ \/^bg- The optimum energy threshold is found to lie between 200 and 700 GeV 
depending on the model ; for our estimates we use two common values, 200 GeV and 500 GeV, 
for all three models to allow direct comparisons between models. 

Figure ^ (left) shows the length of observation time needed for a 3a detection of the FSR 
signal from each dwarf for each annihilation model, for an energy threshold of 200 GeV. Figure ^ 
(right) shows the minimum integrated flux above 200 GeV that can be detected at the 3 a level in 50 
hours of observation for a few ACTs, and the flux predicted by each model for various dwarfs. It is 
clear from these plots that the uncertainties in the astrophysical factors, and therefore observation 
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Figure 2: Left: Observation times needed for 3(7 detection with MAGIC parameters. The three points for 
each dwarf correspond to (from left to right) models A, B and C with a 200 GeV energy threshold. Vertical 
bars correspond to the uncertainties in the astrophysical factors of the dwarfs, with the circles corresponding 
to central values. The solid horizontal line denotes 50 hours of observation time with parameters of the 
MAGIC telescope. The dashed line, at 500 hours, is equivalent to 50 hours of observation with an order of 
magnitude increase in the effective area, i.e. CTA parameters. Right; Integrated fluxes above 200 GeV; the 
dot-dashed, solid, and dashed lines correspond to approximate sensitivities of VERITAS, MAGIC, and CTA 
respectively, for 50 hours of observation time. 



times and sensitivities required for a positive signal, span several orders of magnitude; as a result, 
no model is ruled out, and no dwarf is guaranteed to give an observable signal. These estimates 
only allow us to conclude that, for current and future ACTs, detection of FSR from dark matter 
annihilation from the above dwarfs is likely, but not guaranteed. Further astronomical observations 
of the dwarfs should reduce the uncertainty in the L factors, allowing more precise predictions to 
be made. 

3. Model Identification : Case Studies 

We investigated the prospects of identifying the dark matter model and its parameters based 
on an observed FSR gamma ray signal, and the effects of changes in the energy threshold, the dark 
matter mass, the hadron rejection capabilities of the instrument, and the strength of the signal. 

Our approach for these studies is as follows. Assuming one of the models (A, B, or C) is 
realized in nature, we generated a set of random "data points" distributed in energy in accord with 
the theoretical predictions of this model, incorporating the energy resolution of the instrument. The 
total number of data points corresponds to the prediction of the model for a particular source and 
telescope parameters. These data points were then binned, choosing the bin width to be approx- 
imately double the energy resolution of the ACT being considered. The dark matter mass will 
in practice be unknown, but the final bin can always be made large enough that Ey = can be 
assumed to fall in this bin. 

For each data set, two background samples, corresponding to the ON and OFF regions, were 
generated. The number of background events in each bin was calculated as the difference between 
the event counts in the ON and OFF samples in that bin. A fit to the binned data points was 
performed with all three models, varying the parameters in each model to minimize x^/d.o.f.. The 
three optimal /d.o.f. values obtained in this way were compared to each other, and the model 
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with the smallest value was declared the "best fit" to the data. We performed this procedure for 100 
randomly generated data sets for each of the models A, B, and C. The analysis was performed for 
two sets of ACT parameter values: Agff = 10^ cm^ and 15% energy resolution, representative of the 
reach of current instruments, and Agff = 10^" cm^ and an energy resolution of 10%, representative 
of the reach of future instruments. 

3a or 5a detection — As an example, we provide the results for the case when the FSR flux 
is on the threshold of detection. We performed this analysis assuming the weakest possible signals 
that can be detected at 3 a or 5 a significance above 500 GeV. Table |I] shows the fit results for 
current telescope parameters; results for future telescope parameters are similar. 



"True" model 


model A,B,C 
as best fit 


Best fit WIMP 
mass (GeV) 


Best fit 
cross-section (pb) 




3 a detection 
model A 
model B 
model C 


71, 12, 17 
2, 62, 36 
0,49,51 


3092 ±531 
3041 ± 524 
3051 ±503 


1235 ±265 

1236 ±275 
1238 ±650 


1.17±0.71 
1.74 ±1.62 
1.26 ±0.77 


5 a detection 
model A 
model B 
model C 


83, 9, 8 

2, 64, 34 

3, 17, 80 


3146±433 
3117±471 
3065 ±489 


1087 ±182 
1148 ±238 
1005 ±629 


1.32±0.76 
1.51 ±0.97 
1.44 ±0.95 



Table 1: Fit results for the weakest possible signal detectable at 3(7 and 5(7 levels by the current ACTs. 

The first entry, the distribution of "best fit" models, contains three numbers, corresponding 
to the number of data sets for which models A, B, and C respectively gave the best fit. For the 
other parameters we also list the statistical error bar reflecting the variation of the best-fit values 
and the jd.o.f . of the best fit among the 100 data sets. The fit results show that the overall 
success rate for correctly identifying the model in these scenarios is 61% and 75% for 3a and 5a 
detection respectively. The success rates of correctly identifying the signal as a 2-body or 4-body 
final state are excellent: 90% and 93% respectively. Best fit values obtained for the WIMP mass 
and annihilation cross section are fairly accurate. 

The reader is referred to for detailed results for all cases studied. Lowering the energy 
threshold from 500 GeV to 200 GeV did not seem to have any clear benefits. For a dark matter of 
a lower mass (1 TeV instead of 3 TeV), there was some improvement in distinguishing between 2 
and 4 body annihilation channels for current telescope parameters. Improving the hadron rejection 
capabilities of telescopes also appeared to improve model identification efficiency (from ~ 70% 
to ~ 82%). Overall, the success rate for proper model identification was found to be fairly robust 
with respect to changes in energy threshold, WIMP mass, energy resolution, and hadron rejection 
capabilities of the telescope, and appeared to depend mainly on the signal to background ratio. 

All cases considered in our analysis shared the following common features. (1) The l\i chan- 
nel was very clearly identified because of its 'edge' feature. (2) The two 4-body final state channels 
were not easily distinguishable from each other because their spectra are very similar. (3) Best fit 
values for the dark matter mass and annihilation cross section were in excellent agreement with 
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their "true" values. The mass was better reconstructed with future telescope parameters because of 
superior energy resolution. The accuracy of the mass and cross section determination was particu- 
larly impressive (of order 1%) in the 2/x channel, again presumably due to the sharp edge feature at 
Ey = m),. (4) Future telescope parameters showed a clear improvement over current parameters in 
terms of correct identification of the model. This was mainly due to the increase in effective area, 
which results in a greater number of events and a better significance. 

4. Concluding remarks 

If the PAMELA, Fermi and HESS anomalies have their origin in leptophilic dark matter an- 
nihilation, current and near-future ACTs have an excellent chance of observing the accompanying 
final state radiation from dwarf galaxies. Unfortunately, lack of precise knowledge of the distribu- 
tion of dark matter in the dwarfs makes the signal flux predictions highly uncertain. If a signal is 
observed, the measured gamma ray spectrum can likely be used to identify the correct annihilation 
channel and dark matter mass — a general conclusion that holds for a range of signal strengths, 
dark matter masses, energy thresholds, and instrument parameters — paving the way to a better 
understanding of the microscopic nature of dark matter. 
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